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ABSTRACT 
Demand for ground water and the resulting need to expand our understanding of ground 
water systems has led to the extensive use of models of several different kinds. 
Limitations of physical-type and analog-type models have resulted in the development of 
more powerful numerical techniques. 

The Trescott, Pinder and Larson (1976) model, one of the most versatile 
finite—difference models available, has been modified to run on the University of Alberta 
Amdahl V8 computer. The model has then been used to re-evaluate the field data for an 
extensive pumping test conducted on an aquifer in east-central Alberta. The test, with a 
pumping phase exceeding 34,000 minutes in duration, was thought to be of sufficient 
length to identify boundary and recharge effects. Simulations have been run in order to 
test the validity of several different conceptual models with the intent of selecting one 
most appropiate for the simulation of a twenty-year production period including 
appropiate recharge effects. 

A favorable comparison between real and simulated drawdown in pumping and 
observation wells could not be improved upon to any great extent by the inclusion of 
various forms of recharge. It would seem that the field test does not supply all data 
necessary to effectively quantify the recharge to the ground water system. 

The effectiveness of a pumping well at intercepting recharge, which occurs as 
direct precipitation on an adjacent recharge area, is also evaluated using the Trescott et al 
(1976) model. The influence of transmissivity, storativity, rate of recharge and pumping 
rate have been considered utilizing a hypothetical ground water basin. Stress, in the form 
of pumping, has been added to the ground water system at steady-state and the results 
analysed. Varying parameters such as transmissivity and storativity may result in an 
increase or decrease in the efficiency with which a pumping well can intercept ground 
water moving through an aquifer system, while varying rate of recharge has little effect 


on such efficiency. 
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|. INTRODUCTION 

The increasing demand for ground water and the competition for the resource 
provides an important incentive for expanding our knowledge and understanding of 
ground water systems. One area that continues to receive attention by researchers is 
modeling, where predictive techniques have been used for aquifer evaluation since the 
early 1900's (Thiem, 1906). 

There are a variety of different kinds of models, which may be categorized as 
conceptual, physical, analog or mathematical. Generally the conceptual model is used to 
represent, in a generalized way, our understanding of how various components of a 
ground water system actually function. The formulation of a conceptual model is usually a 
preliminary step in the development of physical, analog or the more powerful 
mathematical models. The physical models are scaled down replicas of real systems. A 
sand-box model for example, may be made up of layers of of sand representing 
aquifers and clay representing confining layers, scaled in proportion to an actual field 
situation. Physical models tend to be very limited in their applicability to practical 
problems but are often very useful for illustrative purposes. Analog models such as the 
resistor —capacitor network type or the conductive-paper type are based on an analogy 
between the flow of ground water and the flow of electric current. These models have 
been used extensively for a variety of practical problems and are particularly useful when 
computer facilities are not available. 

The mathematical approach involves describing, in mathematical terms, response 
of an aquifer system to pumping for example. Governing differential equations can be 
solved using a variety of mathematical techniques to provide a description of drawdown 
or hydraulic head distributions within a region of interest. The analytical techniques, which 
are based on classical mathematical techniques, usually can be applied only to relatively 
simple field situations because they often are based on assumptions such as an infinite 
areal extent for the aquifer, radial flow conditions, and constant aquifer parameters. 

The computer-based procedures, which involve a numerical solution of the flow 
equation, represent a powerful set of tools for solving practical problems. The numerical 
techniques are not subject to the limiting assumptions, which characterize the analytic 


models. For example, they can incorporate a variety of complex boundary conditions, and 


sai Ranntidigernas 98 yarn ciipliy GHabsaney Vorwbein treme VO 12 yigr@y vee 
ax beau Bi leGetn asad BAT YORE Imnrsterertra/” °c > Dom SaaNOR 

« Ko-emrmnoares aochinv INGA Ie pArlonsne war 1 6 eavenees’ 

w hie ar peom antaboros a So aleMMnIt tT ovo? ales ensign 

unos tq Sern eet pele ediewla *s Tors. jeg eae WaT Stee 

A arene nay ho deailos+ WnbeaGo oe “ebormoagte win ont eter, 

gudregn v6) alter is Yo peyal Yo qu pew ad with alas’ %Ot Ineo area 

pre? Wise fa OT nomisqe iy ou Sees FI)" Sin? Oa Qrseewigs* yal 

Hiei eo? MANIA Cee "HOT A aeons yay OF Ob Phat e2bOm 

emiae foe Biebart Qosiws avsouws eyliser Glu. *o) oribee yey ‘aye om 

¥ Vial re cement NS mCIYS etReP Heine. ooo egyt Fewer OtoamEe= 9 

avert sighare kane? dnevus civgela *o wo! orb b” Felewe owen ‘o wot aa 

here hetney gression aie tne ernst 9 tandmeng To) Viney & 35) ibvldnetya a 

Gs stu Pod ae sag 

anon arm? ncametoe nO gridiaeb. covlovns tars tas esirervigr ape aT 

ad ous aholmupn lsaranattib pW Aeyoe? Siorrrxe ior ghia mejaye wtlbae na te’ _ 
nivopwias te foiidnmnet © piiiorh + onvcininet hoitartarttion to" anny s Ore Gee: 
Nouri. ‘ ; att oF keDit el Jeet Lic. somes &ruthy 2 ‘abud agi bger Siuaaet ate 
irgters S2'ylene Peed ein eeupiniagt | trim inate 

MURAG @e Oss Ow NANO Vers eth =e¢ > aig OG 


“se ~ethigdane Ste RATIOS Wel) bet solupE AO TT MOM 
a a 

. joe mT aViGuh) ToIIW 29 W/hagO 1g ee end Bryan _ 7 
Idare Pee Grigg wl 2lod to tea Wawona : "a 


a (ite -roltgnuRed grit! elton eaidia | 
Aen os (Oo VEY & BImeqoaA nes yar? alae 


permit aquifer parameters to vary as a function of both space and time. In addition, it is a 
simple matter to consider multiple discharge or recharge wells. Such flexibility and 
generality is extremely difficult or impossible to attain with either analytic or physical 
models. 

Numerical models are easier to set up than resistor—capacitor models, avoiding 
the selection of appropiate resistors and capacitors and the painstaking soldering of 
numerous connections. Advantages of the numerical approach are evident for 
three-—dimensional situations, which are particularly difficult to simulate with an analog 
model of either type. Changes to a numerical model are simple in comparison to the 
analog or physical model. Such versatility is particularly useful during the model 
calibration phase of the study involving matching of simulation results to the field results. 

The successful application of predictive numerical techniques to real problems 
requires sets of parameters characterizing the system under consideration. While a good 
deal of information is available for the common parameters such as transmissivity and 
storativity, characterization of the quantity of recharge to an aquifer system has been a 
longstanding hydrogeological problem. To effectively utilize the ground water resources 
of an aquifer and to predict the results of various withdrawal schemes, it is necessary to 
have some knowledge of the rate at which the aquifer system is being recharged as well 
as the effectiveness of the particular withdrawl scheme in intercepting recharge. 

Recharge from precipitation varies both spatially and temporally. This variation has 
made the quantification of recharge difficult although a variety of attempts have been 
made to determine how much recharge occurs in a given area. Published rates of 
recharge, expressed as a percentage of the local annual precipitation, vary from about 2 
percent for a small park-land area near Devon, Alberta (Farvolden et al., 1963) to over 
18 percent for a sandy area in Goshen County, Wyoming, U.S.A. (Rapp et al., 1957). Toth 
(1968) estimated recharge near Olds, Alberta to be about 9 percent. 

This study is concerned with computer modeling of aquifer systems. It's goals are 
to adapt the Trescott et al.,(1976) model to the University of Alberta Amdahl computer, 
to re-evaluate a detailed aquifer test conducted in an area near Kirkpatrick Lake, in east 


central Alberta, and to use the model to gain a better understanding of aquifer recharge. 
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The specific objectives of this study with regard to flow modeling are to: 
obtain and interpret available hydrogeological data from the study area and 
formulate a conceptual model of the ground water regime, 
use parameters obtained from conventional interpretative procedures and general 
field observations to construct and calibrate a numerical model by comparing 
computer predicted results to existing field results, 
simulate various pumping strategies with a view to establishing the probable 
long-term response of the system and 
evaluate, with the numerical model, how various ground water withdrawal schemes 
and aquifer parameters influence the interception of water moving through an 
aquifer system when recharge results from direct precipitation on adjacent areas. 
The approach to meeting this objective is based on a sensitivity analysis of a series 


of hypothetical cases. 


Il. GEOLOGY 

The purpose of this section is to present a detailed description of the geology of 
the Kirkpatrick Lake area (Fig. 1). Information presented here has been obtained from 
lithological and geophysical logs available from the oil industry as well as private and 
government publications, reports and maps. Upper Cretaceous strata of the Bearpaw and 
Belly River Formations comprise the near-surface bedrock. These rocks are 
uncomformably overlain by the Pleistocene glacial deposits. A generalized stratigraphic 
column for this area is presented in Figure 2. 

The Bearpaw Formation is the uppermost bedrock unit in the entire study area. It 
consists of an interbedded sequence of marine shale and sandstone. Rocks near the top 
of the sequence crop out south and east of the study area along Sounding Creek and 
Monitor Creek respectively. Test drilling indicates that the Bearpaw Formation is 
approximately 150m thick although the total thickness is not represented in the study 
area because the upper part has been removed by erosion. The exposed Bearpaw strata 
are typically weathered to a light brown color and are very soft. lronstone concretions 
are common (Irish, 1967). 

The Bearpaw Formation has been divided into two members in the study area (Fig. 
2). The upper member is referred to as the Paintearth Member and the lower member is 
referred to as the Young Creek Member (Lines, 1963). The Paintearth Member is 
comprised of sandstone which may be argillaceous, bentonitic or glauconitic. The 
sandstone zones are typically separated by grey shale. At several locations the boundary 
between the two members is a layer of chert pebbles (Lines, 1963). 

According to Williams and Dyer (1930) a definite and traceable sand member 
exists about midway in the Bearpaw Formation. This very fine to medium-—grained sand, 
referred to as the Bulwark sandstone, is variable in thickness generally ranging from 5m 
to 7m. The exact extent of the Bulwark sandstone is unknown, although a similar 
sandstone at roughly equivalent elevations is known to exist about 13km south (Williams 
and Dyer, 1930). 

Subsequent to the work by Williams and Dyer, water well drillers began to report 
that the Bulwark sandstone could produce either white or brown water which suggest 


that it was not a single sandstone but was at least two distinct units. More recent work 
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Figure 2. Generalized stratigraphic column for 
the Kirkpatrick Lake area (Adapted from 
Gliissoldy 1972+ Shaw and Harding, 
1965; Lines, 1963) 
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by Curry (1968) and Clissold (1972) has demonstrated the existence of at least 3 distinct 
sandstone units in the study area (Fig. 3). As a result, the usefulness of "Bulwark" as a 
stratigraphic unit is questionable. Clissold (1972) referred to the upper, middie, and lower 
sandstone units of the Paintearth member in the study area as zone |, zone Il, and zone Ill 
respectively to avoid futher confusing use of the term Bulwark. 

The Young Creek member of the Bearpaw Formation is predominantly grey shale 
or silty shale with some argillaceous or bentonitic sandstone (Lines, 1963). Bentonitic 
beds ranging from a few centimetres to as much as Om in thickness have been reported 
in this member (Given and Wall, 1971). The base of the Bearpaw Formation is defined by 
a layer of chert pebbles (Lines, 1963) and apparently rests uncomformably on 
carbonaceous shale of the underlying Belly River Formation. 

The Belly River Formation in general consists of an interbedded non-marine 
sequence of pale grey sandstone and buff-—colored sandy shale (Irish, 1967). The 
interfingered character of the predominantly non-marine Belly River Formation. and the 
underlying marine Lea Park Formation results in as many as ten members, which are 
arbitrarily placed in the Belly River (Shaw and Harding, 1954). Thus, the Belly River 
Formation in east-central Alberta and particularly in the study area appears to represent a 
mixed marine-continental origin. The continental or deltaic deposits forming the 
predominantly sandstone members such as the Brosseau, the Victoria and the Ribstone 
Creek appear to have been derived from the west or southwest and thin to the east and 
northeast (Shaw and Harding, 1954). Tongues of predominantly marine shale such as the 
Shandro, the Vanesti and the Grizzly Bear members present in the study area generally 
thin to the west. Minor coal seams or coal stringers are present where the Belly River 
Formation becomes more continental in nature. 

The distinct sandstone and shale members of the Belly River Formation, which are 
present in the study area, are shown in Figure 2. Individual members have been described 
by various authors such as Allan, (1918) and Shaw and Harding (1954) as follows: 
Brosseau member — fine-grained, grey calcareous sandstone with some sandy, brownish 

shale particularly to the east. 
Shandro member - dark grey shale containing calcareous and arenaceous concretions. 


Victoria member — fine to medium-—grained grey sandstone and brownish-grey, 
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in study area (Adapted from Clissold, 
1972; Shaw and Harding, 1963) 
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carbonaceous, silty shale with local thin coal seams. 
Vanesti member - grey shale, silty shale, clayey shale and fine sandstone. 
Ribstone Creek member —- soft grey sandstone with some carbonaceous shale and coal. 
Grizzly Bear member - dark blue-grey shale with some ironstone and sandstone 
modules. 
Oldman member - pale grey sandstone and buff-colored sandy shale (Irish, 1967). 

The Oldman member, which has formational status in southern Alberta, represents 
the northeastward pinch-out of a typical Belly River facies and as a result resembles the 
undivided Belly River Formation (Shaw and Harding, 1954). 

Information on the Belly River Formation in the study area must come from 
drillhole information because the unit does not crop out. Lithologic logs from existing 
drillholes and an. electric log from an oil company test hole about 15km to the north and 
on strike with the bedrock of the study area indicate that the Belly River Formation is 
approximately 100m thick in the area. 

Surficial deposits in the general area may be classified as residual deposits, 
unsorted glacial deposits, reworked glacial deposits or transported deposits of alluvial, 
lacustrine or colluvial origin (Wyatt et al., 1938). Very little of the unconsolidated material 
in the study area is residual with most deposits being glacial in origin and subsequently 
reworked by wind and/or water. The Viking moraine (Warren, 1937) occurs as an 
elevated feature north-east of the study area and much of the alluvial gravel, sand, silt 
and clay forming part of the flats around Kirkpatrick Lake has been derived from this 
topographically elevated feature. The Viking moraine is composed of unstratified clay 
with few pebbles or boulders that often comprise a terminal moraine on the prairies of 
Alberta. 

The surficial deposits in the study area are mapped as ground moraine, alluvial 
deposits, aeolian deposits or lacustrine deposits as shown in Figure 4. The ground 
moraine is generally thin with a surface characterized by scattered pebbles and boulders 
of Precambrian origin. Occasional patches of gravelly outwash occur. The alluvial 
deposits tend to be well-sorted with some sand and gravel lenses. Some boulders may 
be present. Both the aeolian and lacustrine deposits are fairly well-sorted with few 


boulders. Bedding may be present in either type of deposit. 
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The sequence of surficial deposits range in thickness from O to 12m (Clissold, 
1972). However, a water well driller’'s log from 16-21-34-9-W4 reports 40m of 
unconsolidated deposits. This thickness suggests that a bedrock channel, perhaps 
preglacial in origin, exists in the eastern portion of the study area. 

The fence diagram shown in Figure 5 presents in detail the spatial relationships 
among the various stratigraphic units. The study area itself is characterized by the 
presence of a broad bedrock high. The bedrock low to the east may limit the lateral 
extent of some of the sandstone units particularly the uppermost sandstone. 

Bedrock units in the study area strike northwest- southeast and dip at a low angle 
to the northeast. This interpretation is based on drill hole data compiled by Clissold (1972) 
for the base of the uppermost sandstone aquifer in the study area The dip direction is 
opposite to that normally expected in this part of the Alberta Syncline (Elliott, 1960) and 
may represent a localized flexure on one limb of the broad syncline. This structure could 


limit the areal extent of the aquifers. 
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Figure 5. Fence diagram constructed from drillhole 
data (For lines of section see Fig. 6) 
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Ill. HYDROGEOLOGY 

In this section, the potential ground water yield and chemistry for various 
geologic units in the study area are described. The data have been obtained from water 
well drillers’ reports on file at the Groundwater Division of the Alberta Research Council, 
from ground water resource evaluation reports prepared by consultants in support of 
ground water licence applications as well as from complaint investigation reports written 
by members of the Groundwater Division of the Alberta Research Council. 

According to Clissold (1972) all test holes drilled in the study area encountered 
aquifers which, individually, or collectively produced ground water in sufficient quantities 
to meet domestic requirements. Yields from wells completed in coarse-grained surficial 
deposits such as outwash sand and gravel range up to 1.0x 10-3m?/s although yields of 
6.0x 10-*m?/s have been reported (Kunkle, 1962). The higher yields likely represent short 
term pumping capibilities. In general, ground water yields from surficial deposits are 
poorer than average where the deposits are thinner and are better than average where 
they are thicker. 

The chemical quality of the ground water from surficial deposits aquifers is 
generally considered poor with total dissolved solids, hardness and sulphate content 
exceeding Canadian drinking water quality guidelines (National Health and Welfare, 1979). 
Total dissolved solids may range up to 2400 mg/I with hardness and sulphates ranging up 
to 600 and 1400 mg/I respectively. The higher total dissolved solid content occurs 
where surficial deposits are thinnest (Borneuf, 1978). 

The Bearpaw and Belly River Formations contain the main bedrock aquifers in the 
study area. Extensive sandstone units up to 20m thick are present in both formations. At 
some locations east of the study area, sandstone in the Belly River represents the only 
reliable bedrock source of ground water. However, the depth of the Belly River 
Formation and the apparent lack of permeable strata in at least the upper part make it of 
minor significance as an aquifer in the study area (Clissold, 1972). Typically, Belly River 
aquifers are found at depths greater than 140m. 

Information on the ground water potential of the three sandstone members of the 
Bearpaw Formation has been available for many years as the result of ground water 


investigations for domestic supplies. The most extensive investigation of the ground 
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water potential of these three sandstone aquifers has been conducted by Clissold (1972). 
The purpose of that work was to locate a long-term supply of water for secondary 
petroleum recovery. The three aquifers are known to be of limited lateral extent but are 
more or less continuous across the entire study area. However, the upper and lower units 
are subdivided by shale units. 

The study by Clissold (1972) covered an area of approximately 200 km? 
immediately to the east of Kirkpatrick Lake. The work verified the existence of the three 
Bearpaw yeni at depths of less than 125m. As previously mentioned, Clissold (1972) 
proposed informal names for each of the units with zone | representing the upper 
sandstone in the area, zone Il representing the middle sandstone and zone Ill representing 
the lower sandstone. Clissold’s (1972) terminology has been adopted in this work to 
identify the various sandstone units in the Bearpaw Formation. 

Zone | was most intensely investigated by Clissold (1972) because early in that 
program there were indications that only zone | could yield the quantity of water required 
for injection purposes. The top of the Zone | sandstone is intersected at depths ranging 
from 10 to 34m. It has a maximum thickness of 20m. The sandstone strikes NNW-SSE 
and dips at about 1m per 5km to the northeast. Cross-sections constructed in east-west 
and northwest-southeast directions show that zone |! subcrops beneath the surficial 
deposits under Kirkpatrick Lake and thins to the northeast and south (see Figures 2 and 6). 

The chemical quality of the ground water from the bedrock varies both laterally 
and vertically. Although no analyses of ground water within the study area are available 
from the Belly River Formation, two such analyses are available for well-water further 
south and east. One analysis is for ground water obtained from the Birch Lake member 
and one is for ground water obtained from the Ribstone Creek member. Analyses 
indicate ground water has relatively high Cl- ion concentration, high alkalinity and a total 
dissolved solids content of about 2000 mg/l. 

The chemical quality of the ground water from zones |, Il, and Ill is shown in Table 
1 along with representative compositions for two aquifers in the Belly River Formation. In 
general bicarbonate, total alkalinity and chloride contents increase with depth, while 
sulphate ion concentrations decrease. It appears that zone Il ground water is better in 


quality than that from either shallower or deeper aquifers. 
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Figure 6. Testholes near Kirkpatrick Lake (Lines 
of section for fence diagrams shown) 


TABLE 1. Ground water chemistry from various 
aquifers in the study area 
(Clrssold, 1972+ *Ktnkle, (962). 
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BR* = Belly River Formation, Ribstone Creek Member 
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IV. EXISTING AQUIFER TEST DATA 

Existing aquifer test data for the Kirkpatrick Lake area consist of results for 
short-term tests by water well drillers, which are carried out in conjunction with 
domestic well construction, and more extensive tests by Currie (1968), LeBreton (1969) 
and Clissold (1971 and 1972) related to evaluation of ground water resources for 
industrial purposes. The most significant testing in the study area was conducted by 
Clissold (1972) in an attempt to locate and evaluate ground water supplies for secondary 
petroleum recovery operations. A total of twenty-nine constant-rate aquifer tests of 60 
to 1499 minutes pumping duration were conducted on zones |, Il and Ill prior to a 
long-term constant-rate aquifer test of 34,510 minutes duration in zone |. 

Fifteen of the 29 constant-rate aquifer tests evaluated the water production 
Capibilities of zone |, 7 evaluated zone Il, 5 evaluated zone Ill and 2 evaluated the potential 
of the Belly River Formation or combined zone Ill and Belly River aquifers. The 15 aquifer 
tests on zone | were distributed throughout the area being studied in order to determine 
the extent of zone | and to obtain information which would allow the selection of the 
most appropriate site for the production well and subsequent long-term pumping test. 
Twenty-six of these tests included water level measurements during recovery. One 
step-drawdown pumping test was also conducted using a well completed in zone |. This 
test was designed to determine the most appropriate pumping rate for the zone | aquifer. 
The results indicated that a properly constructed well could be pumped at rates as high 
as 1.06x10-?m3/s for a short time at least with drawdown being within acceptable limits 
(Clissold, 1972). | 

Graphical analysis of the pumping phase data and, where possible, the recovery 
phase data using the Jacob Method (Cooper and Jacob, 1946) are included in Clissold’s 
(1972) report. (See Tables 2 and 3 for a summary of the relevant data). Transmissivity 
values for zone |, the most extensively tested and promising aquifer ranged from 
2.6x 10-5 to 4.7x 10-3 m?/s for the pumping phase of the test and from 3.8x 10-5 to 
3.4x 10-3 m?/s for the recovery phase. A map showing aerial distribution of transmissivity 
based on interpretation for the early part of the drawdown data and the early part of the 
recovery data is shown in Figure 7. This map indicates that the transmissivity increases in 


a north-easterly direction from Kirkpatrick Lake (Clissold, 1972). A water well driller’s 
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log of a domestic well in 16-21-34-9-4 indicates that the zone | aquifer is thin at that 
location and yields poorly when bail tested. The transmissivity is thus assumed to 
decrease further to the northeast based on the qualitatively low value at that well. 
Because no observation wells were used with any of the short term pumping tests on 
wells completed in zone | no storativity values are available. 

The static water levels in the wells were measured, usually after the test holes 
stood undisturbed over night, and a map of non—pumping water levels has been prepared 
(See Table 2 and 3; Figure 8). This map shows a trend of decreasing water level northeast 
from Kirkpatrick Lake. This flow pattern may be an indication of recharge from the lake 
to the zone | aquifer, which subcrops beneath the lake. 

During the long-term pumping test water levels in 11 of the observation wells in 
addition to the pumped well declined in response to the pumping (Clissold, 1972). Of 
these 11 wells, 10 were completed in the zone | aquifer and the remaining one was 
completed in the zone II aquifer. Although this well was completed in zone Il it was at a 
site where a dual yet separate zone | and zone || aquifer completion was attempted in the 
same well. The fact that the water levels were very close to the same elevation in zone | 
as measured in the outside casing and in zone Il as measured in the inside casing suggests 
communication between the two aquifers, perhaps as leakage along the casing. This 
contention is further supported by water level data from other sites. When separate, 
closely spaced wells were completed in zone | and zone Il the water elevations were 
different and tended to decrease with increasing depth (Clissold, 1972). It therefore 
seems reasonable to assume that the water level decline in the well completed in zone Il 
is the result of a poor casing seal and is not a response to pumping the zone | aquifer. 

Barometric pressure was monitored during the pumping test and small water level 
changes in observation wells were corrected for the influence of atmospheric pressure 
variations during the test (See Table 4). 

Two additional observation wells appeared to respond to long-term pumping 
although no consistent direct relationship was apparent. The water level in the well 
completed in the aquitard above zone | slowly rose prior to the test and continued to do 
so during the entire test. The rate of water level rise varied and it was not until the 


pumping phase of the test was concluded that the water level began to level off and 
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TABLE 4. Observation wells in which water 
levels declined during pumping 


Well number Status yaa Drawdown? 
Saal D2 0 18.24 
SB 74 O* Onal Hore ul 
3C-—i.1 O 6) Sao 
171-71 fe) a) 2h gee) sha TS 

2 Se fe) OS 2047 326 
214 Os fe) O525./ We 2S) 
SA-7 | fe) 1646.3 1.86 
Sy fe) 3256 OV Eras 
3) ey fe) SinzelerS 0.61 
TEZ fe) 4039.6 D208 
1O=7-1 fe) 4436.0 Omir 
F eal al 6) 557983 Orla 


Distance from pumping well (m) 
Drawdown (m) 


Toa 
Drawdown? 


P? = Pumping well 
O* = Observation well 
12-71s* = Shallow well 
12-71d* = Deep well 
TF’ = Thulien farm well 
0.22% = Drawdown corrected for atmospheric pressure 


changes 


" 7. j 
; ae! > A 7 ; 4 
ts ' fb o 
602 e. . 
ee a Y 0 
GS. ¥. ic 
~.  ¢ é 
‘ce .0," ra @ 
face es E 0 
: ; 2] 
\ ee 9 5 
Sr prerer- o 


. . ioe 5 
euzasg sireddsomsaeaot. deen) ioe 
aopigic 


| 


Tihs 


24 


finally began to decline. 

The water level in well 3-71, completed in zone Il, rose during the first few 
minutes of the test then began to slowly decline reaching the pretest level after 600 
minutes. For the remainder of the test, the water level fluctuated within a few 
centimetres of the pretest level. The water level decreased about 5 centimetres within 
20 minutes of the termination of the pumping and then began to rise slowly. 

The analyses of the drawdown data from the long-term aquifer test consisted of 
time-—drawdown plots and evaluations for each well and distance—drawdown for 8 well 
pairs after a pumping duration of 30,000 minutes. A summary of the transmissivity and 
storativity values calculated by Clissold (1972) from the time-drawdown data is 
presented in Table 5. 

The time-drawdown curves indicated that the hydraulic characteristics of zone | 
appear to vary as a function of time. The data summarized in Table 5 are taken from the 
early part of the drawdown curve in each case and thus represent the hydraulic 
characteristics in the immediate vicinity of the test well. Figure 9 is a map showing the 
aerial distribution of transmissivity as determined from the early part of the data from the 
long-term aquifer test. The trend of increasing transmissivity towards the northeast from 
Kirkpatrick Lake is not as apparent as it was in the plot of transmissivity values taken 
from the short-term pumping tests on the individual wells (Figure 7). The differences are 
to be expected as the pumping tests on individual wells give a local transmissivity 
representative of the aquifer in the area immediately adjacent to the test well while, the 
long-term pumping test on the centrally located well and the subsequent drawdown 
analyses with individual observation wells result in an average transmissivity of the 
aquifer between the pumping and observation wells. 

Storativity values determined from the early part of the time-drawdown curves 
of the observation wells range from alow of 1.4x10-4 to a high of 1.6x 10-3, about one 
order of magnitude. The areal distribution of storativity is shown in Figure 10. Storativity, 
in general, increases with increasing distance from the pumping well. This increase of 
about one order of magnitude was initially interpreted as leakage but later computer 


simulations suggest this was not the case. 
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TABLE 5. Summary of transmissivity and 
Storativity values obtained from 


semilog and log-log analysis of 


time-drawdown data. 


Well No. Status to i be Si 
3A-71 p? 0 16.76 
3B=771 O- 6 | 19-00 
Sail fe) SA MT Pas, 15.54 4 
71 O 5a s0 46.64 1 
be Se O S525) Zoe el 
f2=71D% O OS 2 ees 
2) Now O 1648.3 3 Gite 
6-71 O 3475.6 6379.18 16 
Sea O SO2 ks S 14.69 
1O=7 1 O 4436.0 14,69 
TES O 4039.6 95.01 
i 7a e) ONL sks 69.914 
r‘ = Distance from pumping well (m) 
T? = Transmissivity obtained from 
Semilog plot (x10° *m+/s) 
T° = Transmissivity obtained from 
Vog- ode plot ol 0g *myis) 
Sife@= SStorarivity #(x10e*%) 
t2-J15\> = Shallowawell 
12-71d* = Deep well 
TF? = Thulien'’s "farm well 
P* = Pumping well 
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(Adapted from Clissold, 1972) 
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V. THEORETICAL DEVELOPMENT AND IMPLEMENTATION 
The purpose of this section is to describe the basic theory of modeling 
techniques and to explain the features of the Trescott et al (1976) model, including 
options available to the user. The following partial differential equation describes general 


ground water flow in a two-dimensional aquifer system. 


oh ahy = 528 ( 1 
alters =eo tans) 


where T= components of the transmissivity tensor in the coordinate directions (L?t-}) 


ea ae 


h=hydraulic head (L) 
S=storativity (dimensionless) 
W(x, y,t,)=volumetric flux of recharge or withdrawl per unit surface area of 
the aquifer (Lt}) 

Equation (1) has been simplified by assuming that the coordinate axes x and y are 
coincident with the principal components of the transmissivity tensor T,, and Tyy 

Although (1) may be solved by analytically, these simplified solutions are applicable 
only to idealized situations with, for example, uniform aquifer parameters, constant 
pumping rates from a small number of wells and simple boundary conditions. Numerical 
solutions using a computer provide the only means of dealing with more complex 
situations, which more accurately reflect the actual hydraulic setting. 

To use a finite difference method to solve the partial differential equation, the 
aquifer is first divided up into a two-dimensional array of rectangular blocks. It is 
assumed that the properties of the aquifer are uniform within each of the resulting 
blocks. The continuous derivatives in equation (1) are replaced with finite difference 
approximations at the nodes centering on each of the blocks. With N blocks the result is 
N equations with N unknowns, which are the hydraulic head values defined in the center 
of each block. 

For all the simulations presented here, the Trescott, Pinder, and Larson (1976) 
model is utilized. It is the most versatile model for evaluating two-dimensional aquifer 
problems. The model can consider a heterogeneous and anisotropic aquifer with regular 
or irregular boundaries. The aquifer can be confined, unconfined or a combination of 
both. The model allows for multiple discharge and/or recharge wells, evapotranspiration 


and leakage from confining beds in which storage effects may be approximated. 
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Model simulation requires the collection of field data in order to define the region 
and develop the parameters characterizing the aquifer system. It is necessary to have 
sufficient data to provide a good understanding of the boundary conditions with regard 
to both type and location and the character of confinement for the aquifer. Generally, the 
entire suite of geologic and hydrologic data is required to interpret these features of the 
system. The magnitude and distribution of transmissivity and storativity values is best 
determined from aquifer tests, which should be distributed throughout the area of 
interest. The aquifer tests in addition to verifying the aquifer type, should also provide 
estimates of possible leakage from any confining beds . It will be possible to develop a 
conceptual model of the aquifer system at this stage, which presents a unifying picture 
of the parameters of the actual field system. With the addition of a grid,the system can 
be defined in a form necessary for input to the computer model. 

Once the data set runs successfully, its validity can be partially verified by means 
of a history match, that is, a comparison of model predictions with actual field. 
observations of a pumping test. Should the match not be as good as expected or desired 
it may be improved by varying one or more of the parameters such as transmissivity or 
storativity, keeping them within reasonable limits, and re-running the model. With a 
calibrated working model it is possible to make predictions of future aquifer behaviour in 
response to various aquifer stresses. 

The specific steps in the simulation process are as follows: 

if Prepare a map of the aquifer complete with location and type of all 
boundaries. Although the actual boundaries of the aquifer system may lie 
beyond those of the project area, they should, if possible, be included in the 
model. Where actual boundaries are very distant it may be necessary to 
specify imaginary boundaries. 

a Locate all pumping and/or recharge wells on the map of the aquifer. Also 
locate any streams or lakes which may be recharging the aquifer. 

3: Prepare a block-centered grid with any wells located as close to their actual 
positions as possible. They should also be as close to a node as possible. The 
ratio of Ax ./bx : and Ay -/ bY. should not exceed 


j-1 
1.5 to avoid convergence problems. The grid should be oriented such that 


Teas a) souiw (itineacte’ bre 
1 We i.ore TuorQuo tar ch nace rise ered 
deivesn ele Biko re eats seth art gripsiony = =r noitiots. (i alded whe 

@ qaiviat cr sidiczan ed. iwi 2bec peiniteecs (eum mio)? eiauat ate . 
psig Sry TION G srngae Ny reais agate He purave whe ort ie 
nce hereveartt big ¢ To nibh sete i weyers ye wa tet ald NG 
lseorm seer oi A hie vnergenenst ene 
swarr et beitoey yisitwo ee ‘reas nitty Hi Wi Sac guy shay Mee tae or 900 
frei? \atrote tv zholeaitieus isha {OOM HGS & ,D 1% stem na ne 

peusel ic batoeqne £2 Hook £6 9c! TON citer: we hival saat gums tos 

no WIVisayaOe tal Howe etoiwenbian sit ta Bipen 1¢ srs qoivaey yet hevQngran wel yl 
BOI. jeyomn SAY Shine) ane aronit eleien caer Mie mary Ceiegenal 

fm Welvered sehape eum? Sa erkityibeic all ea qitiedind 2: % wba gmitiow 8 

arta us 9 tats BUOSRY OF pmecual 


- - 


ewolor 2s ee ResTONS ORs re wt ni Ces singe oT \ 

“4 4 

ic 10 Gye brensdsoclritiw glaigmoy wiere en co gna T 7 
© 


gil Ygrt Mmetayew Hae. elt 12.29) Who? eoips Prt Cena saree 
ont © Sebulon aq sideese N pee yur? Gire Cage ot 16 ara? oneyed 7 
of vragen sed vere an syaarenis dy Sees Wh aye utes are vena 
- evMehrnbd verge ytigege 
Quis Aptivpes eA lo ger er't no aliawonp miioe> 2c bom gry ly SiGe Ye 
ra Nah cb. qnignir ae 4feum (ow ool oo armcoets yin etesel 
tanto er pio kn atic USN ye Ay bing Beousnan eld em - 
| WSDON'E OF Hae Sh OG OF e Ohiewle yorlt ridisaca an enouinag 
onl 


- aur “AS Sra en dae dines \o citar 
acted ENOKEOG ger aerieitas x sonagieanirs baw ae! 


_ 


30 


there are a minimum number of nodes located outside the aquifer boundaries. 

4. Determine T and S values and their distribution from field data and plot them 
on the grid system superimposed on the map of the aquifer. For an 
anisotropic aquifer the grid should be oriented parallel to the principal 
direction of the transmissivity tensor. 

3), Initial head conditions must be selected prior to model simulation. Often the 
changes in head resulting from a simulation are adequate and the computed 
head values are not important. In such cases the model does not need to start 
from steady-state conditions in which flow is occurring. The model does 
however have the capability of treating more complex problems of head 
distribution. 

S Prepare a data set of the initial hydrogeological parameters to be used for the 
computer model. 

7. Type the data set and run the program. 

Specifically, the Trescott et al (1976) model uses a block-centered, variable size 


grid such as that in Figure 11 and equation (1) can be expressed as the approximation: 
1 Ah ee Ne 1 A dh 
= Trx— <j) Poc—— — Tyy— -( Tyy— 
az ( Brel \ ott) yea "ay, [ 77 Perley 
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Em eA ta 2 
where x = space increment in x direction for column j (L) 


y = space increment in y direction for row i (L) 


i = index in the y direction 


j = index in the x direction 
k = time index 


W(x,y,t) the source term is computed in the model as: 
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Figure 11. Representative block-centered variable-size 
finite difference grid with index scheme 
forenode (iy) From «TrescotGletval., «1976 
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, 
ly ds i 
a (i.5,%) = evapotranspiration flux per unit area (Lt-}) 


flux per unit area from a confining layer (Lt-}) 


A confining layer, stream—bed or lake—bed from which leakage is taking place and 
the effect of storage as being considered is approximated by an equation modified after 


Bredehoeft and Pinder (1970). The equation is as follows: 
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where hi. . = hydraulic head in the aquifer at the start of the pumping period (L) 
»Jx9 
Ze = hydraulic head on the other side of the confining bed (L) 
IV» 
Kr, Z = hydraulic conductivity of the confining bed (Lt"}) 
i Fis: thickness of the confining bed (L) 


s {£,7] = specific storage in the confining layer (L~) 
(Kr gtlm , 8. le i) = dimensionless time 
t = elapsed time of the pumping period (t) 
In this model, evaporation is considered as a linear function of depth below the 


surface of the land by using the following equation: 
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where Ses = maximum evapotranspiration rate (Lt) 
Hn depth below the surface at which evapotranspiration ceases (L) 


G. .= elevation of the land surface (L) 
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Three methods of numerical solution are available for the Trescott et al (1976) 
model. These are: 
1. Line successive overrelaxation (LSOR) 
2. Alternating—direction implicit (ADI) and 
2. Strongly implicit procedure (SIP) which are in order of increasing complexity 
with the more complex numerical solution generally allowing for more rapid 
convergence and the solving of more complex problems. 
All three methods make use of matrix solutions to solve the algebraic equations which 
are generated for each node. 

Trescott et al (1976) concluded that SIP is the more powerful of the three 
interative techniques for most applications. It is particularly useful because it converges 
quickly and does not require the parameter maniputation that the ADI and LSOR 
techniques do to obtain optimum performance of rate of convergence and accuracy. SIP 
does have a trade-off for its apparent benefits in that it requires more arrays and hence 
more computer storage. 

The model used in this study has been adapted for the Amdahl 470V/8 computer 
employing the FORTRAN G compiler and using the SIP method of numerical solution. The 
source code and fixed-dimension arrays require 100K bytes of memory assuming all 
solve routines are compiled with the FORTRAN G Level 21 compiler. This can be reduced 
to 88K bytes by compiling only one solve routine. If all options are included, storage 
requirements are (100 + X/256) K bytes where X is the dimension of vector Y in the 
main program. 


Y dimension = (15 + Na) Nx x Ny 


where Na = total number of arrays required for the options as outlined in Table 6 


Nx = number of columns 


number of rows 


Ny 


Examination of drill hole lithologs and the time-drawdown pumping test data from 
the test holes completed in the zone | aquifer indicate that in the study area the aquifer is 
confined. The extent of the aquifer has been determined from the drill-hole data and by 


deductive reasoning. Drill hole data suggest zone | subcrops beneath the surficial 
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TABLE 6. Number of arrays (Na in equation (6)) required for 
various model options. 


MODEL OPTION NUMBER OF ARRAYS 


Water table 

Conversion (also requires arrays for water table option) 
Leakage 

Evapotranspiration 

SIP 
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deposits in the vicinity of Kirkpatrick Lake in the western part of the study area 
Information obtained from oil company test drilling to the north indicates that the aquifer 
extends well beyond the limits of the study area in that direction. Drill hole 8-71 and the 
Thulien well (10-12-34-9-4) intersect the zone | aquifer so it can be assumed that it 
extends beyond the eastern limit of the study area although how far exactly is not known 
at this time. Zone | tested poorly during the bailing of the well at 16-21-34-9-4 and at 
12-14-33-9-4, so it can be concluded that zone | has thinned to the northeast and to 
the south at these localities and may represent the limits of the aquifer in these 
directions. These assumed hydraulic boundaries occur well within the study area The 
aerial extent of the zone | aquifer within the study area is shown in Figure 6. For 
purposes of the simulations, the western, southern and northeastern boundaries are 
considered to be no flow boundaries. The northern and eastern boundaries are imaginary 
no flow boundaries as the actual boundaries may occur well beyond the limits of the 
study area. However, in most cases these two imaginary boundaries are sufficiently 
distant from the pumping well to prevent them from influencing the model simulation. The 
grid system developed for the area to be modeled is shown in Figure 12. The node 
spacings have been kept small in the vicinity of the wells, gradually increasing away from 
the well yet not exceeding a spacing factor of 1.5 (See Fig. 12) for adjacent nodes as 
suggested in the model documentation. The larger spacing away from the wells reduces 
the number of nodes required and yet does not seriously affect the accuracy nor limit 
the proposed use of the simulation results. 

Transmissivity values determined at the various test holes from the early 
time-drawdown data are plotted on a map of the aquifer, contoured (See Fig.7) and used 
to prepare a grid of transmissivity values for the simulation. Early time-drawdown data 
are used as initial grid values because the transmissivity values assigned to the grid must 
not reflect the influence of any significant recharge or barrier boundaries, or the 
boundaries would have a secondary influence when added explicitly to the model during 
simulations. The transmissivity values, upon which the model simulations were ultimately 
based range about two orders of magnitude, from 0.02x10° to 2.3x10-*m’/s (See Figure 
13). In general, the value of transmissivity increases towards the north and the northeast, 


away from Kirkpatrick Lake. 
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Grid system for the study area (small grid 
widths not numbered) 
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Transmissivity for the model (x10~*m?/s) 


Figure 13. 
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The values used in the preparation of a grid of storativity (See Fig. 10) could only 
be determined at locations where pumping tests were conducted with observation wells. 
The storativity grid used in the model was prepared using only the storativity values 
determined from early results of the long-term pumping test The values ultimately used 
in model simulations which vary approximately one order of magnitude, from 6x10-* to 
80x 10-*, are shown in Figure 14. Storativity, in general, is lower in the vicinity of the 
production well and to the southeast It tends to increase to the southwest, the northwest 
and the northeast 

The transmissivity and storativity values on Figure 13 represent those values 
which were established through trial and error variation in an attempt to obtain the 
parameters which would resuit in the best correspondence between computer predicted 
drawdown and actual drawdown. The values ultimately used for model simulations varied 
less than one order of magnitude from those derrived from the analyses of field data A 
comparison of Figures 7 and 13, and Figures 10 and 14 indicate that in general, 
storativity values had to be decreased and transmissivity values had to be increased in 


order to obtain the best match between field and model! data 
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Model storativity (x10~‘) 


Figure 14, 
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VI. MODEL SIMULATION AND RESULTS 


A. Kirkpatrick Lake Study Area 

This section describes first, the comparison of analytically and digitally derived 
solutions, second the various models, which are each potentially representative of the 
hydrogeological conditions in the study area, and third the application of the model to the 
investigation of the efficiency with which a well can intercept water moving through an 
aquifer. 

Before any attempts were made to utilize the simulation model, tests were carried 
out to ensure that it had been modified to operate on the University of Alberta Amdahl 
computer. The test data used to verify the model are shown in Figure 15. Boundaries for 
the domain are assumed to be impermeable, and transmissivity and storativity values are 
constant for the non-leaky, artesian, test aquifer. The boundaries are located sufficiently 
far from the pumping well that the aquifer can be assumed to be infinite for the time of 
simulation. The model verification is carried out by comparing results of the simulated 
drawdown at two observation points with the drawdown calculated from the analytic 
solution for the same problem. Early—time divergence between analytically derived 
drawdown values and values determined by computer simulation is greater than late—time 
divergence. This early-time divergence is consistent with results obtained by others 
during similar simulations. Prickett and Longquist (197 1) suggest that this error results 
from discretizing time for successive time steps during simulation. Thus, calculated 
drawdowns are somewhat distorted during early parts of the pumping phase. Accuracy 
can be improved if necessary by decreasing the size of the time increment. However, in 
practice it is convenient to place less emphasis on results from at least the first few time 
steps. The close correspondence of the calculated drawdown from the two methods 
indicate that the model is operating satisfactorily. (See Table 7, Figures 16 and 17). 

Attempts to select field conditions on which to base simulations suggest, initially 
at least, that five different but equally possible arrangements may be appropriate. As a 
result, five sets of conditions are simulated during modeling in an attempt to evaluate the 
possible responses of the system to pumping and provide a basis for estimating the 


most likely long-term behavior of the system. The five scenarios for recharge can be 
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TABLE 7. Comparison of analytically calculated 
drawdown with computed drawdown 
ee UG W(u)? re Calculated Computed 
drawdown(m) drawdown(m) 
yoo 1402 Sal 76.8 es is! P00 
1.41 56 4.61 1.46 1.40 
rede Be ois fale, 1.08 17.0 
Ney fe, o.00 woo 1.2.80 
7.42 11 6.24 36 PANO E 
Us (2) Giead, Bawes 2.14 222 
1351.0 4.3 Ueeulst 2.20 Fe eN 
2s U0 220 160 Zee 2.40 
OS bere OOS) 425 12.2 Oi0n 0.00 
1540 S5 0.24 0.08 ne 
2s 64 50 O56 Omen? O20 
4.57 743) 0.98 0.30 O230 
7.42 18 leer 0.42 0.40 
bis63 Wa 1.74 Oro 5 tr elere 
13210 7.4 26 10 Vio Oa) 
27.19 4.8 2a 0.80 0.80 
t' = Time since pumping began (min.) 
u? = Argument 
W(u)? = Theis well function 
r* = Distance from pumping well (m) 
4050 = "x10 * 
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summarized follows: 


Model 1 - Withdrawal of ground water from a system where recharge is negligible. 

Model 2 - Withdrawal of ground water from a system where induced recharge from a 
lake is approximated by a constant boundary. 

Model 3 - Withdrawal of ground water from a system with leakage from the overlying 
confining shale unit. 

Model 4 — Withdrawal of ground water from a system with the lake a constant head, and 
leakage from the overlying confining shale unit. 

Model 5 — Withdrawal of ground water from a system where recharge from the lake is 
controlled by the hydraulic conductivity of the sediments of the lake bottom. 

The five model systems differ from each other only in the type of recharge 
conditions. Recharge to the aquifer can be modeled either using boundary conditions 
providing flux or leakage through confining units. 

Model 1 represents a case where no recharge is assumed to occur. The south, 
west and northeast boundaries are considered to be impermeable and hence hydraulic 
barriers. Those to the north and east are located sufficiently far from the pumping well 
to effectively remove their influence during simulation. 

Model 2 represents a case where recharge from Kirkpatrick Lake is assumed to 
occur. The main difference between model 1 and model 2 is that in model 2 the lake is 
considered to be penetrating the aquifer where it subcrops and 's thus in direct hydraulic 
connection with the aquifer. During simulation the lake, which is the western boundary, is 
represented as a constant head boundary and may recharge the aquifer along a 5 km 
front, once the cone of depression reaches the boundary. The southern and northeastern 
boundaries are considered to be impermeable, and the northern and eastern boundaries 
are located sufficiently far to effectively remove their influence during the simulation. 

Model 3 represents a case where recharge in the form of leakage is assumed to 
occur. The difference between model 3 and model 1 is the leakage from the confining 
shale unit which occurs in model 3. The confining shale unit has been assigned a vertical 
hydraulic conductivity of 10-4! m/s which could be typical for a bentonitic and silty shale 
(Kruseman and DeRidder, 1970) such as occur in the study area. Model 3 and model 1 


have identical boundary conditions. 
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Model 4 represents a case where recharge from the lake, which is represented 
as a constant head as well as from an overlying confining shale unit, is assumed to occur. 
The boundary conditions for model 4 are similar to model 2. Additional recharge is now 
possible as leakage from the overlying confining shale. Model 4 is essentially a 
combination of models 1, 2 and 3. 

Model 5 represents a case where leakage can occur as induced infiltration from 
Kirkpatrick Lake with actual quantities controlled by the hydraulic character of the lake 
bed sediment. The aquifer in which the pumping well is compieted subcrops beneath the 
lake and, for the model 5 simulation is considered to be in hydraulic connection with the 
lake. Field investigation indicates that the mucky, organic rich lake bottom sediment is 
predominantly silt and clay. Initially, because of lack of actual field data, the lake bottom 
sediment is assigned a thickness of one metre and a hydraulic conductivity of 1x10-’ m/s 
for simulation purposes. The southern and northeastern boundaries are considered to be 
impermeable and therefore hydraulic barriers. Those to the north and east are located 
sufficiently far to effectively remove their influence during simulation. 

The calibration of the model is intended to demonstrate which set of hydraulic 
parameters and boundary conditions satisfactorily represents the actual hydrogeologic 
conditions which occur in the field. The five sets of model systems are designed to 
represent a spectrum of possible situations that should include the actual conditions. 
Because of the good quality of data available for the area, it is felt that most of the 
uncertainty in representing the systems rests with the recharge conditions. 

The computed drawdowns for each of the five models are compared (See Table 
8) to the actual drawdown values obtained from a 24 day pumping test, conducted at a 
rate of 1.06x10-?m3/s, during which the water levels of 11 observation wells, completed 
in the same aquifer as the pumping well, were monitored (Clissold, 1972). Figure 18 
presents, in graphical form, the variation between observed and mode! 1 computed 
drawdown values. Except for one obvious point, corresponding to observation well 
3C-71 located about 6 m. from the pumping well, the match between simulated 
drawdown values and observed drawdown values is satisfactory for each of models 1 to 
5. These results represent the best match after attempts made to improve the match 


between observed and computed values by trial and error varying of transmissivity, 
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TABLE 8. Summary of actual drawdown and drawdown 
computed by models 1 to 5 for the 
pumping and observation wells after 24 
days of pumping 


Well No. DD' Computed drawdown (m) 
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storativity and leakage within reasonable limits. 

The simulated drawdown at the pumping well exceeds actual drawdown by about 
1.5m (8%) for a simulation in which recharge is negligible (model 1). In simulations which 
have some recharge (models 2,3,4 and 5), the simulated drawdown at the pumping well 
exceeds actual drawdown by about 0.6m (3%). The actual drawdown at the closest 
observation well (3B-7 1) exceeds the computed drawdown by about 0.1m, or less than 
1%. Simulated drawdown for observation well 3C-71 is consistently much higher than 
the actual drawdown for all five cases, and as a result some problem with the original 
field data is suspected. Actual drawdown values for observation wells 11-71, 12-71, 
9-71, 6-71 and 8-71 deviate about 10% from computed values in all five cases. The 
general trend of the computed drawdown values is similar to that of the actual 
drawdown for observation wells 10-71, 7-71 and Thulien’s farm well located 
approximately five kilometres from the pumping well. The results do however, deviate up 
to 0.1m from the actual values, an error approaching 50%. Attempts to improve the 
match were unsuccessful. No model seems to represent the actual hydrogeologic 
conditions any better than the others and as a result does not yield computed drawdown 
values which fit the observed values any better, at least not after a period of 24 days of 
pumping. 

Original field data (Clissold, 1972) are no longer available to review in an attempt 
to explain the anomalous drawdown value for observation well 3C-7 1. It may be 
significant that when the observation wells were installed, well 3C—71 was the only one 
with which a preliminary short-term pumping test was not conducted. Although the 
reason for the lack of a test test could not be confirmed, it is possible that some 
difficulties were encountered during the construction and development of the well. One 
possible explanation for the apparent anomalous behavior of observation well 3C-7 1 
could be related to the fact that it has been completed at a location where the bedrock is 
the shallowest observed in any of the test holes drilled. The shallow bedrock could have 
been fractured by the over-riding glacial ice of the Wisconsin ice advance resulting in an 
area where recharge from the surface occurs very readily. Very rapid recharge to the 
well through a fracture network could result in less drawdown in the observation well 


than the simulation predicts. 


eae FRET MAD Nedeyane 
eartt Geel 10 comiimaaree ei art atuanme (fT 
iit arigietaunn (rewiterahas Bh PE -ENaW haters SY wren 
\snight> art divw mstddie ering Huser a ex bins 20-85 aut De Were 
‘vast 17-1 ( ahew nonaioade i) Zvliv-worawe tS feuiA Seragene a ete 
enT <teea dvit ev bevy Bopemes lest 0! ructa relives 1 TB itne 188 
surat ort! *S tat or ralrmsa) 2elayirrurwe bec tt to kent Ma 
aTe30) isi ye eS rts [Tee LOT shane eileen sot mwobwanb 
dw Bibel” tevayrorody aides ait lew gageriga orl mags savamals evlt yierOIROIGA 
bt svvane of SqmaHA AOE grsa}eiqan sare ve eoutty arse ee nowt mr Oot 
Sigh lanpr thy levis ert hh ie dmeae ighom olf utreegmumry sow tame _ 
wee web onhgieae biely ion aa 10Z%) & ‘al Sig ert: afl cwety “untied: yaw seven 
, ae 6S te boivweaq « yelieton tensl?s carted viw esuiky Osu auto one Bt tite: nadae? 
| sarigeriagy 
samate hs PUNBIVaL OL aldpieva Sagnal Co ene (STET Ateasto) emer elaM larign® 
| ed-yerten. [SDE law neakwmEdo TO° gully Cowedwnty auc. welt signe OF 
ioyiric orth ww 1 -DEMew nelieti@ ay) allen norevwRed at marty teat mn 
aut NUON beToyBnos tarlasW. teat prligrhug mrret—rora yarirenianny & sairtwy AW 
Bmon tut sigiesog, 21) Hernmyinas ad "69 blecs sen) tee? 6 te Wome of) vo? nossa), _ 
At) low oft Yc InerngolevekBos HotsunyAgs Yt HONS vOm=ingoOTe ew aBiUOD: 
(5235. Hew) Aoitavisedo to tolvatied eucinmans mormon att ReREaGRD sidesog 
® doorbad sr ead rictisas! ss DeTalQmos riged ont mrt toed ert oF barele7 ed barge 
west Yoorboe wane onl bein eaiont worett te yr a Laviaadto tamwotirta 
| — ani menooalW St Te.65) qoniy reve arty Reo ond, | 
eget ia pay Broo koe nwe ane mot agranrgsn overt ean 


. 


50 


Figure 19 shows a comparison of observed and computed drawdown using 
model 1 after 4, 8, 12, 16, 20 and 24 days for the pumping well and observation wells. 
It should be noted that the correspondence between observed and computed drawdown 
after pumping of the production well for periods of 4, 8, 12, 16 and 20 days is about 
the same as that for the full pumping period of 24 days. Similar results, although not 
included here, are achieved using models 2, 3, 4 and 5. 

The results of the simulation trials show that each of models 1 to 5 can equally 
well simulate the hydrogeologic response of the aquifer system at least up to a time of 
24 days. It appears that the aquifer test was not of sufficient duration, when combined 
with the sufficiently low recharge rate, to establish whether Kirkpatrick Lake will interact 
with the aquifer system over the longer term. Yet, a 24 day pumping test would be 
considered by most hydrogeologists to be a relatively long-term test. These results call 
attention to the general problem of being unable, with a single test, to evaluate pertinent 
boundary effects. 

In order to evaluate the long-term response of the aquifer to pumping at a rate 
of 1.06x 10-*m3/s, each of the five models was used as the basis for a 10 year 
simulation. Table 9 is a summary of the decline of water level in the pumping well with 
time for each model. Using model 1 of the system (no recharge), calculated drawdown 
exceeds the available head of 21.4m in the pumping well in just one year of pumping and 
reached a total decline of 42.68m in ten years. With model 2 there is a reduced rate of 
total water level decline, although drawdown also exceeds that available in just one year 
with a total decline of 26.23m after ten years. Model 3 exceeded available drawdown in 
one year, ultimately declining to a level of 38.83m, not significantly different from results 
of model 1. With model 4, drawdown in the pumping well exceeds available drawdown in 
one year but resulted in the least total drawdown of any of the models used, 25.75m at 
the end of 10 years simulation. With model 5, drawdown exceeds that available in one 
year and gives rise to the second lowest total drawdown after a 10 year simulation, 
27.1m. Model 5 is selected for a longer-term simulation in order to determine the time 
required for steady state conditions to develop. Steady-state conditions are attained by 
year 16 with no further changes in head at later times. The maximum drawdown was 


27.1m resulting in a water level which was 5.7m below the upper limit of the aquifer. 
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It is unfortunate that the calibration technique did not verify that one of the five 
models considered during simulation would be more appropriate than the others for use 
in predicting the effects of continued ground water withdrawal and hence be suitable as 
a ground water resource management tool. Model 5 was selected for the 20 year 
simulation as this model, qualitatively at least, is possibly the most appropriate for the 
following reasons. Firstly: the static water levels in test holes completed in the zone | 
aquifer declined away from Kirkpatrick Lake towards the north-east suggesting that 
recharge from the lake is a distinct possibility. Secondly, the study area is characterized 
by saline soils indicating considerable evapotranspiration and if the Thornthwaite method 
is used, the area is characterized by having a net deficit of soil moisture (LeBreton, 1969) 
and, theoretically at least, little or no recharge. Hence, model 5, with leakage from 
Kirkpatrick Lake the only source of recharge to the zone | aquifer could reasonably be 
selected for long-term simulations of the effects of ground water withdrawal in the 
study area The 20 year time period has been chosen for the extended simulation as this 
was the expected duration of the pumping of ground water for injection purposes. Long 
term production of ground water at the rate proposed will, based on model predictions, 
result in the decline of the available head to the point where dewatering of the aquifer 
will take place. At this time it is likely that the rate of decline of the water level will 
accelerate as the effective transmissivity decreases. This decline may be slowed or 
halted once the cone of depression has increased enough in size to begin accepting 
recharge from Kirkpatrick Lake. However, by this time the cone of depression will be 
large enough to interfere with the water supplies of domestic users in the area The 
severity of interference can only be evaluated when the available heads for the wells 
being intefered with are known. Thus, it seems likely that long-term production of 


ground water at the proposed rate of 1.06x10-*m’/s is, to say the least, impractical. 


B. Theoretical Analysis of Recharge 

Long-term ground water management strategies require some understanding of 
the effects various aquifer parameters have on the efficiency with which a well (or wells) 
can intercept recharge. The remainder of the simulations in this section are designed to 


evaluate how effective a pumping well can be in appropriating recharge in an active 
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ground water flow system. The basic approach involves a sensitivity analysis of the 
response of a model of a large flow system to pumping. With recharge explicity defined, 
it is possible to relate it, with mass balance calculations, to the source of water 
withdrawn from the aquifer. Figure 20 describes the hypothetical ground water system, 
which forms the basis for this group of simulations. The grid system is comprised of 
twenty rows and twenty-five columns with a uniform distance between nodes of 1000 
m. The five zones shown on Figure 20 represent five areas with differing rates of 
recharge. The actual magnitude of the recharge can be varied as desired for steceesve 
simulations. Zone | is part of the basin where no recharge is assumed to occur. Zone 3 
has a recharge rate which is double that of zone 2. Zone 4 and Zone 5 have recharge 
rates which are three and four times respectively that of zone 2. 

The hypothetical basin has no-flow boundaries to the north and west 
representative of a ground water divide and constant-head boundaries to the south and 
east, representative of a stream at those locations. Recharge in the form of direct 
precipitation is applied to the zones shown in Figure 20. As expected for steady state 
conditions, the hydraulic gradient declines away from the no-flow boundary in the zone 
of greatest recharge towards the constant—head boundary representing the stream. 
Stressing an aquifer in the basin by inserting a pumping well should result in the 
formation of a cone of depression which may intercept ground water moving through 
the aquifer possibly resulting in steady state head distribution. Comparison of the amount 
of water withdrawn from the system by the pumping well during the simulation with the 
amount of recharge introduced to the system as the result of direct precipitation, once 
steady state has been reached, should give an indication of how effective a well can be in 
the interception of recharge for that particular set of hydraulic conditions. The more 
effective the pumping well is, the greater will be the percentage of water derived from 
recharge. 

All computer simulations involve the Trescott et al., (1976) model. For this type of 
simulation, it is necessary to input the hydraulic head distribution in the natural flow 
system as an initial condition. Thus, computed drawdowns are superimposed on the 
natural flow system. The first step in this simulation procedure is to compute initial 


hydraulic head distribution by setting storativity to zero, omitting any stress, such as 
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Figure 20. Hypothetical ground water basin 
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pumping, on the aquifer and allowing the simulation to proceed to steady state. Once the 
pre-pumping head distribution is obtained, it can be used as the initial head in the 
transient simulations comprising the second step in the modeling procedure. 

The hydraulic parameters required for the transient simulations include recharge 
rate, storativity, transmissivity and pumping rate. The set of preliminary values, selected 
arbitrarily for the first simulation, are designed to be representative of actual geological 
systems. Transmissivity is 1.14x10-*m?/s, storativity is 1.0x 10-5, and recharge has been 
set at a maximum (that is, zone 5 on Figure 20) of 4.0x10-"m/s, equivalent to ten percent 
of an assumed annual precipitation of 50 cm. Pumping rates are allowed to vary from 
8.33x 10-4*m3/s to 133.0x 10-*m3/s. Simulations using the initial hydraulic parameters 
resulted in steady state conditions in a few days or a very few weeks at the most, even 
with the higher pumping rate. To obtain data from which conclusions might be drawn, the 
initial hydraulic parameters are varied in trial and error procedure in an attempt to choose 
those which would increase the time to reach steady state. Ultimately seven sets of 
hydraulic parameters are selected for simulations. Each of the seven simulations is run 
using five different pumping rates in an attempt to determine if pumping rate affects the 
effectiveness of the well at intercepting recharge for a given set of hydraulic 
parameters. 

Tables 10, 11 and 12 present a summary of the simulation results, describing 
percent of recharge intercepted by a hypothetical pumping well for five different 
pumping rates and seven different sets of hydraulic parameters. 

Storativity values are increased one order of magnitude for simulations one to 
two with other hydraulic parameters identical. The result is an increase in the length of 
time required for the simulation to reach steady state (See Table 13). The results are 
consistent for each of the five different pumping rates used for these simulations. 

Simulations three and four are another pair of cases that vary only in the value of 
storativity. Again, the simulations with the higher storativity value require a longer period 
of time to attain steady state (See Table 13). The influence of storativity on the 
effectiveness of the well at intercepting recharge is probably explained by the shape of 
the cone of depression. Decreased storativity results in a cone of depression which 


increases in both radius and height. The resulting increase in area over which water 
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TABLE 11. Summary of computed percentage of 
pumped water which is intercepted recharge 
for various pumping rates and hydraulic 
parameters. 
Time Test 4 Test 5 
(days) PR' PR? PR? PR* PR PR! PR? PR? PR* PRS 
1.24 0 0 0 0 0 0 0 0 0 0 
3209 0 0 0 0 0 0) 0 0 0 0 
9407 0 0 0 2 0 0 0 0 2 2 
LO 0 0) 6 0 2 0 0 6 3 6 5 
163.3 0 6 iS 2) 4 13 ee 12 13 T3 
2537 0 Le 2) 9 9 25 24 24 24 24 
3956 t3 18 15 18 17 38 35 36 30 3a 
60.8 25 aS) zy 30 aS) 50 53 52 50 51 
Die eid 38 4] 42 42 42 63 65 64 66 65 
140 50 =) ahs) =i) 07 LD a 76 78 2a 
214 63 qs 70 69 68 88 88 88 87 87 
318 75 82 82 81 Sir 100 94 94 94 94 
479 88 88 88 90 90 100 Si) 97 98 
I19 2) 94 94 96 96 100 100 u8 
1080 100 100 98 98 he 100 
1620 100 she s}e) 
2430 100 100 
3650 
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TABLE 12. Summary of computed percentage 
of pumped water which is intercepted 
recharge for various pumping rates and 
hydraulic parameters. 


1.24 0 0 0 0 0 0 0 0 0 0 
o.09 0 0 0 0 0 0 0 0 0 0 
5.07 0 0 0 0 0 0 0 0 0 0 
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16.5 0 6 3 S 4 0 0 0 2 1 
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TABLE 13. Summary of computed time for a hypothetical 
well to reach steady state for various 
pumping rates and hydraulic parameters. 
Number of years to reach steady 
Simulation state 
No. buy R? Se 08235" (6.06. 2.2 OO. Og meson sce 
1 Jibs ecm eats O40) 3 0.04 O2O7- 2007 rae 
2 aie: PC COS 2 0 25 0.38 O.58— 10.58 0.80 
3 1.14 Live LOe Aaied 2:5 On 38 OS 8eg0 258 0.80 
4 1.14 Ose 225,96 Ze 4.44 6.66 6.66 
S) 220 Sew TO? Pe 0t8'7 Uses AIT HOS oK 20 
6 toe Pe eee 6 Fare Ns) 4.44 6.66 6.66 
a) epee 2 tO ‘6,66 0,00 210.00 NR° NR> 
Tiez2Transmissivity. (xXi1007m-7S) 
R? = Recharge rate (percent) in zone 5 
S2 = Storativity 
O8.co° = Pumping rate (x10 ‘*m77s) 
NR° = Steady state not reached 
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moving through an aquifer system can be intercepted makes it possible for the well to 
intercept recharge more effectively. 

The quantity of recharge was also examined as a variable influencing the behavior 
of the aquifer. In simulations two and three, rates of recharge in zone 5 are ten percent 
and four percent of the annual rate respectively, with identical values of transmissivity 
and storativity. Although initial head distributions are different for each of the two 
simulations, there is no difference in the length of time required to reach steady state 
conditions for a specific pumping rate. Similar results are obtained for cases four and six 
where recharge rates in zone 5 are four percent and two percent respectively (See Table 
13). It appears that varying the rate of recharge, while maintaining other parameters 
unchanged, does not effect the efficiency with which a pumping well can intercept water 
moving through an aquifer system. It is likely that there is a lower limit below which the 
recharge rate cannot be reduced for this result to hold. 

Simulations four and five are based on identical rates of recharge and storativity 
but transmissivities are 1.14x 10-3m?/s and 2.28x 10-m?/s respectively. The simulation 
with the higher transmissivity attained steady state more quickly than the simulation with a 
lower transmissivity for each of the five pumping rates. Similar results are obtained with 
cases six and seven which again differ only in transmissivity, (1.14x10-3m?/s and 
0.57x 10-3m?/s respectively). Consistent with previous results, the simulation with the 
higher transmissivity attained steady state more quickly for each pumping rate. 

Decreasing the transmissivity value while maintaining the same pumping rate 
results in a cone of depression which is increased in height but decreased in radius. 
Decreased transmissivity will also result in a reduction of the rate at which ground water 
will move through the aquifer system. The net result in the simulation when transmissivity 
is halved is a decrease in the efficiency with which a well can intercept recharge. 

Results from simulations suggest that the efficiency of a pumping well to 
intercept recharge is not related to pumping rate when hydraulic parameters remain 
unchanged during simulations. Figures 21 and 22 present the results of a mass balance 
from two of the seven simulations, and indicate that generally, the efficiency of the 
pumping well at intercepting recharge does not vary more than five percent for pumping 


rates ranging over 1500% (8.33x 10-m?/s and 133.33x 10-*m’/s). However, for a specific 
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pumping rate, varying hydraulic parameters results in variation in the length of time for 
the simulation to reach steady state conditions (See Figure 23). One exception noted is 
that if only the rate of recharge is varied (Tests 2 and 3, 4 and 6), there in no difference 
in the length of time required for the simulation to reach steady state conditions. 

The seven sets of test data and the five pumping rates indicate that the efficiency 
with which a pumping well intercepts recharge varies according to the nature of the 
hydraulic parameters. Some variation of the parameters results in increased efficiency, 
some in decreased efficiency, and one parameter, the rate of recharge, at least appears 
to have little effect on the efficiency of the pumping well to intercept water flowing 


through an aquifer system. 
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Vil. CONCLUSIONS 
1. It has been possible to adapt the Trescott et al (1976) finite difference model for 
aquifer simulation in two dimensions to the University of Alberta Amdahl computer and to 
use the model to re-evaluate a detailed aquifer test in the Kirkpatrick Lake study area 
located in east central Alberta. 
2. The model which was calibrated with the data from a 24 day pumping test did not 
effectively differentiate between five different arrangements of hydraulic parameters. 
Although the pumping test utilizing 11 observation wells was of a time duration that 
locally at least would be considered very long, it did not adequately assess the boundary 
conditions nor the recharge conditions. It is apparent that for a particular class of 
problem, especially where the impact of pumping extends over a large region, long 
pumping tests on a single well combined with a large number of observation wells will 
not necessarily yield the best results or give valid information on which to base model 
development, whether analytical or numerical. It would be more beneficial to strategically 
locate a number of pumping wells each requiring at least one observation well in order to 
better evaluate aquifer and recharge parameters as well as to identify and locate 
boundaries. 
3. Since the field test results do not allow proper assessment of recharge conditions, it 
is only possible to use the calibrated model qualitatively for any long-term simulations. 
With the static water elevations in the test holes completed in the zone | aquifer 
decreasing with increasing distance northeast from Kirkpatrick Lake, and the zone | 
aquifer appearing to subcrop beneath the lake, the simulation incorporating leakage from 
the lake may prove to be the most appropriate model with which to predict the response 
of the aquifer to the proposed rate of pumping. A 20-year simulation resulted in a 
drawdown of 27.1 m in the pumping well, reaching steady state at about 16 years. The 
drawdown at the pumping well exceeds the available head by about 25 percent. 
4. Assuming recharge occurring as direct precipitation exceeds some basic minimum, 
varying the rate of recharge does not effect the efficiency with which a pumping well 
can intercept ground water moving through an aquifer system. Decreasing storativity 
improves the efficiency while halving transmissivity reduces the efficiency with which a 


pumping well can intercept recharge when other parameters are unchanged. 


oa" me edad Aa waGgsy ot eee a behetoe & 


smodia 

toa bib deer onigriug yar 6S snow aie eer itive somohas aaw Aorta % 
piacere ailugrey > aieineRa ria Tree YiS ou) renuoried snacenattib 
ine reine ae & foaen Glew neievcecco Tt priite rae gyildjrrasey ett 
Watruod @0 efeeds vistupebsven bib af one! ev D2 °eHiewbe od bluow tant te 
‘> e281) “alUnU dS To} say INE Real etpitieno> Sg marae ont YON 
“eo Eee) $08) @ VC FONSI Ae Brien 7c 'GQhern gf “artw ylaoegee y vale 
ive Hew MOlTee 1eOh To Waly awh ily. Lemndnns i Ses 1 sie@? 
REA Sanh or tigitive Ne HOTem gin bey 4euig tc Wieck) "hd ort? olery vineReesen: 
Wades rus Of iGalteted oie éc mloswe 4! Jase oo Gacylan wee Je 
=? VSbAS ola motevaseds ene Inge! th Teigies rome wiv’ prdatrneg “oe were 2 


ateool ons. yiiieebi. 0) &e Now 26 Cale nadia cpryetor) hrw whine awl 


ry varieties aap e290) Yo WteiTahSEas TeQo7C wale lot ob aruaen Teed iphe) ery won? £ 


ROcHELarmte Pres PaGi ye IOF ype wiBL lanord Ssteubhes ot emeot sidered yo 
Wiis» | eto eto bagionon dec” 7207 ah wh efetiinvans law wee ert? 


brio ee bee MOREE Mee cebarMon amuse Proce tw rere 
reg ee eT SMEG a SIh nl inibecte it, go let on) frare? Gow of pMoags 


- 


eenoqen ys 7H ibang. 27 nai Avy inbooi) 696 Toe seg wet Of oe eyo NT VE Mee 


2 >) ceupery “italy ng ‘Bay te he te eal 416.416). hetent> iq “i a caries arte He a 
> 


on? avtey S1 ods 16 wind: 5 re goiter: vy Oruariwe ae co Ay ‘ce te | or 


+ 
Say | ba 
Prurinn Ks aioe 2oe0a"e NGM NIA “2esib ah pute spuartan: grumyeeA | 7 
lew. wg ra a mw date > Chi Mm afoo aniancey te ete anit pF wv 


rasaiag, 25 ‘vod vy Dole oldaleds. at! thacoss lw qeenpriiey ert te 


67 


REFERENCES 


Allan, J. A. (1918): Sections along the North Saskatchewan River and the Red Deer and 
South Saskatchewan Rivers, between Third and Fifth Meridians; Geological Survey 
of Canada, Summary Report for 1917, Part C, p 9-13 

Borneuf, D.M. (1979): Hydrogeology of the Oyen area, Alberta Research Council report 
78-2, 35 pages. 

Bredehoeft, J.D. and Pinder, GF., (1970): Digital analysis of areal flow in multiaquifer 
ground water systems: A quasi three—dimensional model; Water Resources 
Research, Volume 6, No.3, p 883-888. 

Currie, D.V. (1968): Chevron Standard 13-1-35-10-W4, Oil Field Injection Source Well, 
33 pages. 

Department of the Environment (1974): Hydrogeology of the Hamilton Lake area, Alberta, 
39 pages. 

Elliott, RH.J. (1960): Subsurface correlation of the Edmonton Formation; Journal of the 
Alberta Society of Petroleum Geologists, Volume 8, No. 11, p 324-338 

Farvolden, R.N., Meneley, W.A., LeBreton, E.G., Lennox, D.H. , and Meyboom, P. (1963): 
Early contributions to ground water hydrology of Alberta; Research Council of 
Alberta, Bulletin 12. p 98-111 

F. W. Schwartz Consultants Ltd. (1980): Modeling of groundwater systems at Stony Plain, 
March, 1980. 72 pages. 

Given, M.M. and Wall, J.H. (1971): Microfauna from the Upper Cretaceous Bearpaw 
Formation of south-central Alberta: Bulletin Canadian Petroleum Geology, Volume 
19, No. 2 (June 1971), p 504-546. 

Gravenor, C.P. and Bayrock, L.A. (1955): Glacial geology, Coronation district, Alberta 
Research Council Alberta Preliminary report 55-1, 38 pages. 

Green, R. (1972): Geological map of Alberta; Research Council of Alberta map 35 (scale 
ale26 7,000) 

Hydrogeological Consultants Ltd. (1972): 1971 analysis of groundwater data, Chevron et 
al Project No. 1, Hamilton Lake, Alberta, 32 pages. 


(1972): Groundwater Report, Kirkpatrick Lake area, 


prepared for Mapco. 153 pages. 


(iW tale? neal cart iO Se al -Y Saas owed ese): 
| . | 7 bs a 

qi ccnA ears aie) Hosmer ott SO yolaeBenBen ATEC arate iver at to wenmrmged 
st hoinuel nor RotnoRaa ede Te falinerns eye rues AO3et) La a cras 

HER ==Se 9) of) erm’ Dispos meuoinee ic yreged are 

{ROT A eget bref 14D sonna!s 22 a hanes JN weitesM WA 

,  toilgnincs Witanehr aril te (aoe aane bovere a enptadnenas oe 

('T+68.q ai rivi oak 
MR yrieitte epetyts wiaerniég SRNegd? (OPet nt serra tuto W Wie 
ogee STORET torah - 7 

wanes 4 otter 2) 73941! wt hare ahustp0iM ‘)TS)) A) TW bem AM ar) 
sriyioty yRo08) (mast TEAnriGling T etallue Syne here etyod be nadernet 
. Sbi-0G8 q '7 2! atk! £ OB 
ered py wait reubgnicts veSiAee bicker neti yA sonal bn ag Terevae 

2ege0 SE 1 GF tee raniriiie emects. iigesod een % 
Wiens) BE Quer ainiRlA: $0 "= hN69 Aataded |kre 9° Yo que anignioae NEN F mere 
(OGS.5 HS. 14 
tm Ne anh 40 Stet bid 2 whearot amgelonpantyl 
ehacté edad nodimetf ae. = 7 
we iste) | 


68 


Irish, E.J.W. (1967): Geology of the Oyen sheet; Geological Survey of Canada Map 21-66 
(scale 1 inch to 4 miles) 

Kruseman, G.P. and DeRidder, N.A. (1970): Analysis and evaluation of pumping test data, 
International Institute for Land Reclamation and Inprovement, Bulletin 11, 200 pages 

Kunkle, GR. (1962): Reconnaissance groundwater survey of the Oyen map area, Alberta; 
Research Council of Alberta Report 62-3, 25 pages. 

LeBreton, E.G. (1969): Groundwater availability in the Hamilton Lake area, Alberta: 
Research Council of Alberta Internal Report, 18 pages. 

Lines, F.G. (1963): Stratigraphy of the Bearpaw Formation of Southern Alberta; Bulletin of 
Canadian Petroleum Geology, Volune 11 No. 3, p 212-226. 

Meyboom, P. (1967): Estimate of ground water recharge on the Prairies; Water 
Resources of Canada, University of Toronto Press. p 

Minister of National Health and Welfare, (1979): Guidelines for Canadian drinking water 
standards 1978, 79 pages. 

Nicols, R.A.H., and Wymann, J.J. (1969): Interdigitation versus arbitrary cutoff: Resolution 
of an Upper Cretaceous stratigraphic problem, western Saskatchewan; American 
Association of Petroleum Geologists Bulletin, Volume 53 p 1880-1893. 

Prickett, T.A. and Lonnquist, GC. (1971): Selected digital computer techniques for ground 
water resource evaluation: Illinois State Water Survey Bulletin 55, 62 pages. 

Rapp, J.R., Visher, F.N., and Lettleton, R.T. (1957): Geology and ground water resources of 
Goshen County Wyoming, United States Geological Survey, Water Supply Paper 
1377, 142 pages. 

Thiem. G. (1906): Hydrologische Methoden: J.M. Gephart, Leipzig. 56 pages. 

Toth, J. (1968): A hydrogeological study of the Three Hills area, Alberta Research Council 
of Alberta Bulletin 24, 117 pages. 

Shaw, E.H., and Harding, S.R. (1954): Lea Park and Belly River Formation of East Central 
Alberta: Western Canadian Sedimentary Basin Symposium; American Association of 
Petroleum Geologists, Tulsa, Oklahoma, p 297-308. 

Trescott, P.C., Pinder, GF., and Larson, S.P. (1976): Finite-difference model for aquifer 


simulation in two dimensions with results of numerical experiments; United States 


Geological Survey, 116 pages. 


aha nee eins Saal aman 

| aligue § 

te seh aed ome To AS amt ic yr a) 
Ass! oq £: les grime: 

ha earls eine we wee VSoW tong ta amet ah ) 9 mace 

a sant aniateT eonaninh! sent 16S ‘ nA 


faye gS Finbar at coniiak (UenPew ow veh a © Wrein 
cegag UF Ove aerate on 


Aafuloce® Fails ye Viow ee evans ware ERE k'. foshayey gra ht ae 
qacineen®:itniearacs eee crue Salen Site cava ths ala aR TS: ae 
Seer GME GES oreutcl/ aut coeqaoiest mpiowe? to nade — 
was 1) 10" Reuse ening inipih qurgrees(? Tt, Qe ‘ayaa AT sean 
aermrc Se’ 42 insellge ene vay wre got, “olevlber onones sate 
‘qunrwrday shee Grup ofa ypoiaet (ont) DA oataites orw 2 wake JS, 
RT Vighue WARW .aVVe crolans PENS beret) orien Tres a 
sagen Rat tet 
eons BE Cithlss Metiaa Mim~nétontarstseqete ae wiGer D meitt 
ruts? Congas ar wdie, 50 en veut" at iT en | cs (Geet 4 
sais Fi BS “dial acacia Yo 
eareD ted 1G Hore avi ye DME eT ee GRE!» MR genet Rem A wot 
Sono AOa A AROMRTA OLutauriye Ninee amainomtibes "elberemc: rie area 
YE TSF G-aimonalvO sehit wetenyenicemy ra 
vertiant: xt I6BGen Soren SAINT) SZ rok oe a er 
patess Reteyl serge vans ier (3) iw eno Ses 


ae 


69 


Wall, JH, Sweet, AR. and Hills, L.V. (197 1): Paleoecology of the Bearpaw and Contiguous 
Upper Cretaceous Formations in the C.P.0.G. Strathmore well Southern Alberta; 
Bulletin Canadian Petroleum Geology, Volume 19, p 691-702. 

Warren, P.S. (1937): The significance of the Viking Moraine; Transactions of the Royal 
Canadian Institute , Volume 21, No. 46, p 301-305. 

Williams, M.Y. and Dyer, W.S. (1930): Geology of southern Alberta and Southwestern 
Saskatchewan; Geological Survey of Canada, Memoir 163, p35-43. 

Wyatt, F.A., Newton, J.D., Bowser, W.E. and Odynsky, W. (1938): Soil survey of Sullivan 
Lake sheet; University of Alberta Bulletin 31, 102 pages. (Appendix by John A. Allan 


- The relation of the geology to the soils in the Sullivan Lake Sheet). 


- a err 4 : 7 y re" 
| | te Se 
ee Pee 3} St = 
7 hill (ie na’ 7 ous 
Pia rH zngitmast ari pee aware 
ne oe 
hei ny ee 
agile Fe Vey Mos SEEN) OVW eaety mew 
not, & nil. per aBOQGAy maipaey ROT VE re ween? 


flivae Bre €7 


ie 


fi 
WE of 
4 


